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On ■ Abstract 

^h ; 

"q \ The influence of different chemical potential for different flavors on color 

superconductivity is analyzed. It is found that there is a first order transition 

iq ' as the asymmetry grows. This transition proceeds through the formation of 

bubbles of low density, flavor asymmetric normal phase inside a high density, 
superconducting phase with a gap larger than the one found in the symmetric 

C"-- case. For small fixed asymmetries the system is normal at low densities and 

superconducting only above some critical density. For larger asymmetries 

<^ ' the two massless quarks system stays in the mixed state for arbitrarily high 

densities. 
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A lot of progress has been made recently in mapping the phase diagram of QCD in the 
low temperature, high baryon density region. The main idea advanced by this research is 
the almost unavoidable formation of quark-quark condensates (pairing) and the spontaneous 
breaking (or restoration) of color-flavor symmetries of QCD |[l]||. These ideas may have an 
impact on the phenomenology of neutron stars and of heavy ion collisions at SPS energies 
where compressed strong interacting matter is believed to exist. 

Since lattice calculations at finite density are unfeasible at this date two alternative 
approaches are used to study this problem. One, valid at asymptotically high densities, is 
the use of perturbation theory. The gaps due to perturbative QCD were believed to be 
small 0, exponentially suppressed as in weak coupling BCS theory. It was recently realized 
y-§|, however, that the main pairing interaction in this regime is due to the exchange of 
magnetic gluons. These interactions make the normal phase even more unstable than short 
range forces do and result on a parametrically much larger gap. Extrapolating this down 
to more accessible densities one could expect gaps up to 100 MeV at chemical potentials 
fj, ~ 1000 MeV. It remains to be seen at which density the sub-leading correction (in powers 
of the coupling) become large and invalidate the perturbative expansion. Another approach 
to the problem is the use of phenomenological models such as NJL-like models and the 
instanton liquid model [|TU|-|r7|] . In particular, a large class of four-fermion interaction models 
was considered. Fortunately, for semi- quantitative purposes, they tend to give very similar 
results. This is due to the fact that the formation of the condensate depends only on the most 
attractive channel between quarks at the Fermi surface. Different four-fermion interactions 
that are both attractive in the scalar-isoscalar color 3 channel and of roughly the same 
strength produce similar BCS instabilities. Renormalization group analysis support this 
conclusion and help in classifying all four-fermion operators according to their importance 
to the formation of the gap ]T8"| , |TT| . 

With the help of these phenomenological models the most likely scenarios of symmetry 
breaking were identified []I],^|,|TD|,|T3[ . In the case of 2 massless quarks chiral symmetry is 
likely restored at some critical density and color SU C (3) is partially broken down to SU C (2) 
by the formation of the condensate (2SC) 

<«2?«S> = (sSeg> ~ ^ 3 e<% (l) 

where the a,/?,... are color indices, a,b,... are flavor indices, i,j,--- are spinor indices and 
qL t n are left (right) handed spinors. The condensate in eq. ([[]) singles out the third color 
(green from now on), only blue and red quarks pair up. 

In the three massless quarks case the symmetry breaking pattern is more subtle pO| , |2T . 
The preferred condensate is the so-caled color-flavor locked (CFL) state 

(Qu^ql]) = -(qm^qRj) ~ ^(^5^5^ + K 2 5 ab 6 0a ) ( 2 ) 

The condensate eq. (|2|) breaks SU C (3) x SUl(3) x SUr(3) down to SU c+ l+r(3), leaving the 
ground state invariant under a combined color-flavor rotation. The remaining symmetry 
SU c+ l+r{3) is the same as the one of the vacuum after spontaneous breaking of chiral 
symmetry (SUv{3)). This suggests that there is smooth connection between the hadronic 
phase to the CFL phase. The spectrum of excitations is identical in both cases: the gluon 
octet is mapped onto the vector meson octet, the quarks onto the baryon octet and the 8 
Goldstone bosons onto the pseudo-scalar mesons. 



The real world lies between the 2 and 3 massless quark case since the mass m s of the 
strange quark is neither much larger nor much smaller than the QCD scale. At arbitrarily 
high densities the presence of the strange quark mass is irrelevant but not at lower densities. 
It is difficult then to determine whether the CFL phase extends all the way to lower densities 
and connects to the hadronic phase or whether there is a window in density with the 2SC 
phase. Another difference between the discussion sketched above and the potential real 
world application is the difference in densities among the different quark flavors. In heavy 
ion collisions there is an excess of ~ 15% in the number of down quarks (minus down anti- 
quarks) relative to the up quarks, and zero net strangeness at all. There are also significant 
asymmetries in neutron star cores. 

Here we will analyze the two massless quarks phase in a model that conserves quark 
flavors 

We choose our interaction to be a contact four fermion coupling with spin, flavor and 
color structure abstracted from the one gluons exchange projected on the color 3 channel 
(the 6 is repulsive). Our starting point is then the Lagrangian 

£ - fi a Af a = ^({icr^ + fi)6 ab + 5Tl b )5 a p%l)p b - g^V^^^o-iV^^^ ~ <U<W> ( 3 ) 

where /i = {jj, u + /Ai)/2 is the chemical potential averaged over flavors (equal to one third 
of the baryon chemical potential), 5 = (/i u — /id)/2 is their difference, M a is the operator 
that counts the total amount of flavor a (quarks minus antiquarks) and the spinor indices 
were omitted. The interaction in eq. (|3]) does not mix left and right handed quarks and 
a sum over handedness is implicit throughout the paper. The matrices a^ are defined as 
o^ = (1, —a) for the left handed quarks and o^ = (1, (?) for the right handed ones. The mean 
field approximation, although not arising as a systematic expansion in a small parameter is 
usually reliable to determine the phase structure of the theory. It is exact in the limit of large 
number of flavors (keeping g 2 Np constant) and on the weak coupling limit. To perform the 
mean field approximation we rewrite the model defined by eq. (Q) by introducing a dummy 
field A aa/3b 

£ + fx a J\f a = iplaiiicr^df, + fi)5 ab + 6^)5^^ + —A aapb A* fbSa e a ^ 5 

+ ^aa^lplbAybSaeap-yS + ^\ b ^*s a A aa p b e a p 1&} (4) 

where e a/ 3 7 ,5 = 5 ai 5^ — 5^5^. It is convenient to use the Nambu-Gorkov formalism and 
write eq. (f|) as 

£ + n a Af a = ( ^ i) T o 2 ) M I ^ J + —A aa p b A* bSa e a p jS , (5) 

with 

^-y $* V ^ - n + 6t 3 ) ' [b) 

where <r M = (1,?) for the left handed quarks and <r M = (1, —a) for the right handed ones. 
The matrix r 3 acts on flavor space (and is diagonal in color and spin) and $ is the operator 



&aaf3b — ^-yaSb ^Safi-, (7) 

acting on color-flavor space. 

The mean field approximation amounts to estimating the path integral by its saddle 
value point. Ignoring then the fluctuations on the Aq^j, field we are left with a Gaussian 
integral over the fermion fields that produces the familiar det(.M) = Tr log(.M) term. The 
potential f2(A,/z, 5) defined by 

e -in(A,n,8) = fx>^V^e^ c+IXaMa (8) 



is the expectation value of the Hamiltonian Ti — fi a M a in the state that minimizes 7i — /i a A/" a 
in the subspace of states with a given value of A. In our case 

% f d 4 p 2 

Q(A, /1, 5) = -J TT^tr logM + —A aapb A* yb5a e a p lS , (9) 

where the trace is over spin, color and flavor spaces. The value of A Qa/ % in the ground is 
obtained by minimizing Q(A, fi, 5) at fixed fi. We make an ansatz, corresponding to the 2SC 
state, for the form of the order parameter 

A Qa/ % = -^-2e ab e Q/ 33. (10) 

As mentioned before, this choice breaks color S c (3) down to S c (2). By using a basis in spin 
space that diagonalizes p. a and by doing the similarity transformation 

M^UMU-\ u =[\l 2 ) (11) 

we are left with the task of computing the determinant of 

where 

<5 a/3 = -iAe a/33 (13) 

acts only in color space and e ± (p) = ±p — fi. Using the relations 

tr log ( £ B D ) = Tr hg(-BC + BDB^A), (14) 



$ 2 = -P, (15) 

where P is the projector on the first two colors subspace, and 

Tr flavor log(a + br 3 ) = log(a + b) + log(a - b) (16) 

Tr color log(o + bP) = (N c - 1) log(a + b) + log(a), (17) 



we arrive at 

Q(A, n,S) = $ + £+,- |J (0[2 log((po + 5) 2 - E± 2 (p)) + log((p + 8) 2 - e ±2 (p)) 

+2\og((p - 8) 2 - E ±2 (p)) + log((p - 8f - e ±2 (p))] 
+C, (18) 



where E ± (p) = J A 2 + e ±2 (p) is the energy of the quasi-particle excitations. Performing 
the po integration and fixing the constant C by demanding Q(A, fx, 8) has the correct (free 
theory) value at A = we arrive at 

n(A,^5) = f-y 1 ^E + AnEHp)-p-o(s-E ± (p))(E ± (p)-s)) 

+ \e ± (p)\-p-0(S-\e ± (p)\)(\e ± (p)\-5)} 

+3p + 6{5 + e-(p)){5 + e-(p)). (19) 

The first line corresponds to the contribution of the red and blue quarks, the second to the 
unpaired green quarks and the last one to anti-particles that can be created when 8 < fi. 
The mean particle number density of each flavor is obtained by 

d 
n u + n d = -— fi(A, /i, 8) 

-sgn(e ± (p)) + 6(8 - ^ (p)\)sgn(e ± (p))] + 9(6 + e"(p)) 

d 
n u -n d = — — fi(A,/x,5) 
06 

All three dimensional integrals are cutoff at momenta A. For illustration, we choose the 
values A = 900 MeV, g 2 /A 2 = 4.4. They produce a gap of around A ~ 100 MeV at n = 300 
MeV and 6 = 0. In fig. (1) we have Q(A,[/,,8) as a function of A for different values of /i 
and 8. They all correspond to states with the same value of n u + rid, but increasing values of 
the asymmetry n u — n d - The minimum of the first curve, top to bottom, corresponds to the 
flavor symmetric, superconducting state. As mentioned above, the superconducting state 
consists of paired blue and red quarks of both flavors, with green quarks left unpaired. 

There are two ways of looking at flavor asymmetric matter: one may fix different densities 
or different chemical potentials. Let us first consider the first possibility and imagine in- 
creasing the asymmetry n u — nd from its the initial zero value while keeping the total density 
n u + n d fixed. States with small n u — rid values can be created by adding to this symmetric 
state up-quasi-particles and down- quasi- holes. In the grand canonical formalism used here 
this is accomplished by increasing the value of 8 from zero to some small value. Since there 
is a gap in the spectrum of blue and red quasi-particle, no quasi-particle with these colors 
are created for 8 < A. The flavor asymmetry is made up only with green, up-quasi-particles 
and green, down quasi-holes. If all colors were paired we would see no change in Q at all 
in the superconducting phase until the split between chemical potentials 8 equaled the gap 
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FIG. 1. Potential U in units o/10 10 MeV 4 as a function of the gap A in MeV. By the order the 
curves meet the vertical axis, from top to bottom: i) [i = 300 MeV, 6 = 0, ii) /i = 300 MeV, 5 = 40 
MeV, Hi) fi = 300 MeV, 6 = 77.5 MeV, {n u + n d ) normal = 3.1 10 6 MeV 3 ,{n u + n d ) supeT = 3.91 10 6 
MeV 3 , (n u -n d ) normal = 2.34 10 6 MeV 3 , (n u -n d ) super = 0.90 10 6 MeV 3 , iv) // = 310 MeV, 5 = 83 
MeV, (n u + n d ) normal = 4.31 10 6 MeV 3 , (n u + n d ) super = 3.45 10 6 MeV 3 , (n u -n d ) normal = 1.04 10 6 
MeV 3 , {n u -n d ) super = 2.69 10 6 MeV 3 , v) y. = 322.5 MeV, 5 = 89 MeF, {n u + n d ) normal = Am 10 6 
MeV 3 , {n u + n d ) super = 3.91 10 6 MeV 3 , (n u - n d ) norma i = 1.22 10 6 MeV 3 , (n u -n d ) sune ,r = 3.14 10 6 
MeV 3 , vi) fj. = 310 MeV, 5 = 120 MeV 
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A. The potential can be more easily lowered in the unstable state with A = since there 
quarks of all colors can be created, even with small 5. As a result Q(A = 0, /i, 5) decreases, 
with increasing 5, faster than Q(A ^ 0, //, 5), as can be seen in the second curve in fig. (1). 
At some value of 5 the superconducting and the normal phase are both stable. The value 
of n u — n d is much higher on the normal phase than in the superconducting phase, by the 
reasons explained above. One could imagine that, as n u — n d increases, the system proceeds 
towards a mixed phase of A = and A/0 phases, with the up-rich A = phase occupying 
more and more of the space until the superconducting phase completely disappears. The 
presence of two conserved charges makes things a little more involved. The reason is that 
the total density n u + n d in the superconducting phase is much higher than the total density 
in the normal phase. Thus, even hough by choosing the right amounts of the two phases one 
can arrange so that the overall asymmetry n u — n d be any number between the asymmetries 
in the two phases, the overall total density is always smaller than the initial total density. 
In order to achieve higher values of n u — n d and keep n u + n d constant we need to increase 
also the value of the average chemical potential \i (the result in shown on fig. (1), fourth 
curve to meet the vertical axis, from top to bottom). The value of \x is chosen in such a way 
as to satisfy 



Tl u ■+" Tl d X[Tl u -r 1^d)\super ~r \i- %)\1T"u, "T Tld)\normal 
'•'« '"d •Ey'l'u <^d)\super "i~ V ■£)\P>u <^d)\normali 



(21) 



for some x, < x < 1. 5 is determined by the condition that Q(A,fi,5) should have 
two degenerate minima as a function of A (so phases can coexist), so we are left with 
two equations eq. (|21~|) and two independent parameters (/x and x) so equations eq. (|2T|) 
can, generically, be (uniquely) satisfied. The transition, then, proceeds through a series of 
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FIG. 2. Phase boundaries as a function of the total density (in units of 10 6 MeV 3 . 



mixed states, each one with different values of \x and S (arranged in such a way as to keep 
n u + n<i constant while increasing n u — rid) ■ Each one of these mixed states is formed by 
a slightly asymmetric, high density superconducting with a larger gap than the one before 
the transition started and a highly asymmetric, low density normal phase. The mixed phase 
exists until its normal component occupies the whole space and the total density of the 
normal phase equals the initial value of n u + rid (fifth curve to meet the vertical axis in fig. 
(1), top to bottom). Increasing n u — rid even further is accomplished by increasing 5 again, 
at roughly fixed ii. At this point the superconducting state is not even metastable anymore. 

The different phases as a function of n u + rid and (n M — rid)/{ri u + rid) are sketched on fig. 
(2). The horizontal axis, corresponding to the symmetric case, is in the superconducting 
phase for any density, since the BCS instability exists for any density. The vertical axis 
correspond to the completely asymmetric case, with ii u = —fid- There are only up-quarks 
and down- ant i- quarks present, no pairing is possible so those states are in the normal phase. 
In between these two phases there is a large mixed phase region that occupy most of the 
ri 'U) n d > region. The precise location of the phase boundaries are somewhat sensitive 
to the parameters of the model. The example shown in fig. (2) was obtained with the 
parameters discussed above. The distance between the superconductor-mixed phase is on 
the large side among those obtained with other reasonable parameters. As fig. (2) shows, for 
asymmetries larger than about 20% the system remains in the mixed phase for arbitrarily 
high densities. The phase boundaries in the low density region, where the model considered 
here has no resemblance to real QCD, continue all the way to the origin. Unfortunately 
it is very hard to determine the properties and precise location of the mixed phase with a 
reasonable degree of certainty. The bubbles of the normal phase are charged and, if formed 
in heavy ion collisions, may lead to charge separation that could be observable. 

The color-flavor locked state, expected in the 3 massless flavor case, is expected to behave 
under asymmetries in a similar manner. Actually, since color is completely broken and all 



quark colors have a gap, one could expect the mechanisms discussed here to be even more 
effective than in the 2SC state. In the case of massive strange quarks though, it remains to 
be seen whether existence of "gapless superconductors" pointed out in |E2| may hinder this 



conclusion. A discussion of this question, as well as a more complete treatment of the two 
and three flavors asymmetric matter will be presented elsewhere |23[ . 
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